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ABSTRACT 

We present a measurement of the expansion and brightening of Gl.9+0.3, 
the youngest Galactic supernova remnant, comparing Chandra X-ray images ob- 
tained in 2007 and 2009. A simple uniform expansion model describes the data 
well, giving an expansion rate of 0.642 ± 0.049% yr~^, and a flux increase of 
1.7 ± 1.0% yr~^. Without deceleration, the remnant age would then be 156 ± 11 
yr, consistent with earlier results. Since deceleration must have occurred, this 
age is an upper limit; we estimate an age of about 110 yr, or an explosion date of 
about 1900. The flux increase is comparable to reported increases at radio wave- 
lengths. Gl. 9-1-0. 3 is the only Galactic supernova remnant increasing in flux, 
with implications for the physics of electron acceleration in shock waves. 

Subject headings: ISM: individual objects (Gl. 9-1-0. 3) — ISM: supernova rem- 
nants — X-rays: ISM 



1. Introduction 

The small Galactic supernova remnant (SNR) Gl. 9-1-0. 3 was determined to have an age 
of less than 140 ±30 yr based on a comparison of a radio image from 1985 and an X-ray image 
from 2007 (Reynolds et al. 2008; Paper I). This makes G1.9-f0.3 the remnant of the most 
recent known Galactic supernova. Its X-ray spectrum is highly absorbed {N^ ~ 5 x 10^^ 
cm~^; Reynolds et al. 2009 [Paper III]) and dominated by synchrotron emission with a very 
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high roUoff energy, hu^ow ~ 2.2 keV. Faint thermal emission can be isolated from the N rim, 
with line emission from Si, S, Ar, Ca, and Fe, as well as a 4.1 keV line from ^^Sc produced 
in the decay chain from '^^Ti (Borkowski et al. 2010; Paper IV). 

The mean expansion of Gl.9+0.3 between 1985 and 2007 was estimated to be 16 ± 3%, 
or (0.73 ± 0.14)% yr~^ (Paper I). A comparison of the 1985 radio image to a new 2008 
vlaQ radio image confirmed this result, giving a mean expansion rate of 0.65% yr ^ (Green 
et al. 2008; Paper II). At an assumed distance of 8.5 kpc, consistent with the high X-ray 
absorption, the implied shock speed is about 14,000 km s~^, comparable to that inferred 
from measured line widths (Paper IV). The high velocity, and the symmetric nonthermal 
X-ray morphology, suggest a Type la origin, though this is far from certain. 

We observed Gl.9+0.3 with Chandra again in 2009 for a total of 237 ks, and now have 
the opportunity of examining the expansion between two X-ray observations. Since the radio 
and X-ray morphologies differ somewhat (Paper I), it is important to determine an expansion 
from X-ray data alone, which should give the best indication of the blast wave speed. 

2. Observations 

Our 2007 observation with Chandra using the ACIS-S CCD camera (S3 chip) took 
place on 2007 February 10 for 24 ks, and March 4 for 26 ks. We checked aspect correction 
and created new level 1 event files appropriate for VFAINT mode. No flares occurred during 
the observation. CTI correction was applied and calibration was performed using CALDB 
version 4.3.1. Finally, the data sets were merged and weighted response files created. About 
8000 counts were obtained. We reobserved on 2009 July 13 and July 27 for a total of 237 ks. 
The same procedures were followed. About 40,000 counts were obtained. We verified that 
the data sets were well aligned for both observations by using the fitting method described 
below to compare pairs of sub-observations of the same epoch. We found negligible shifts 
(far smaller than a pixel). The weighted average interval between observations was 2.408 yr. 

We used an elliptical mask, chosen by inspection, to exclude non-source parts of the 
image, allowing the analysis to proceed faster. The mask also excluded circular regions 
surrounding several point sources. We checked that the mask did not skew the results by 
repeating the analysis without it; the results did not change significantly. Therefore, all re- 
sults below were obtained with the mask. In addition, we subtracted a constant background. 
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estimated from a region of the 2009 image free of discrete sources. 



3. Method 

To measure the expansion, we compared the 2009 platelet-smoothed (Willett 2007) 
image to the raw 2007 data. (We tested for dependence on smoothing method by repeating 
the analysis with an image smoothed using the spectro-spatial method of Krishnamurthy, 
Raginsky, & Willett (2010). The results were not significantly different, either in best- 
fit values or in uncertainties.) The model contains four parameters: a physical scaling 
factor, a surface brightness scaling factor, and expansion center coordinates. We shrank the 
background-subtracted image by an overall multiplicative factor about a point chosen by 
eye to be near the center, scaled its brightness (counts/shrunken pixel) by another factor, 
added back a background taken from the 2007 image, and shifted the image as necessary 
to minimize our fit statistic described below. The shift from the originally chosen origin 
produced an "expansion center," just the point on the 2009 image from which this simple 
uniform-expansion model yiel ds the closest approximation to the 2007 data. This method is 



analogous to the one used by I VinkI (120081 ) for Kepler's supernova remnant 



Because the number of counts per bin is small, the 2007 data are not well-represented 
by a Gaussian distribution, making the use of the statistic inappropriate. The data are 
better described by a Poisson distribution: the probability of getting a particular set of data 
counts {Di\ is given by 

M 

L = X{-^eM-M.). (1) 

where = Si + Bi is the sum of source and background model exp ectation va lues. The 
maximum- likelihood fit can be found by minimizing the Cash statistic ( CashlE979l ) given by 



C = 2 5^ [M, - A + A (log A - log M,)] (2) 

i 

In the limit of large numbers of counts, the C-statistic becomes identically x^- The minimum 
was found using a downhill simplex algorithm. The resulting values for the parameters were 
found to be independent of the choice of initial values. 

The method produces a factor that represents a change in the surface brightness (counts/pixel) 
(after scaling by the ratio of exposure times, a factor of 4.779). However, the effective change 
in pixel size as the 2009 model is shrunk to match the 2007 data means that the X-ray flux 
changes by an additional factor, just the linear scaling factor squared. Our best-flt value for 
the linear scale change is (1.55 ± 0.12)%, or an expansion rate of (0.642 ± 0.049)% yr^^. Our 
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brightness scaling factor by which the 2009 model had to be diminished, before correcting 
for the exposure-time difference, was 0.2064 ± 0.0042, indicating a surface-brightness change 
of a factor of 1.01 ± 0.02 over the 2.408-year interval. After correction for the size increase, 
this yields a mean flux increase of a factor of (4.55 ± 2.08)% or (1.89 ± 0.86)% yr~^. The 
derived expansion center is shown as the cross on the X-ray image of Fig. 1. Our results are 
summarized in Table 1. All errors are 90% confidence intervals. 

We also perf ormed Markov ch ain Monte Carlo (MCMC) simulations using the PyMC 



software package (jPatil et al.ll2010l ). These simulations converged on the same optimal values 
of the parameters found using our optimization routine presented in Table 1, to an accuracy 
much better than our reported errors, supporting the validity of the maximum-likelihood 
method. 

By inspection, the remnant appears to have spatial variations in its relative expansion 
factor and change in brightness. The east and west X-ray bright lobes appeared to expand 
by a larger factor than the radio-bright northern region and the dimmer southern region. 
To see if this visual impression was supported statistically, we divided the remnant into 
four regions and performed a joint fit allowing the regions to have separate expansion and 
brightness factors, but sharing a common expansion center. The resulting fit was not a 
significant improvement on the four-parameter fit. Confirmation of spatial variations in the 
expansion rate will require a deeper observation. 



4. Flux Increase 

The flux increase can be measured directly by performing a joint analysis of spatially- 
integrated spectra. We used a simple absorbed power-law model, but instead of subtract- 
ing background we modeled it as a combination of sky and particle components (Paper 
III). Background spectra were extracted from an annulus surrounding the remnant (Pa- 
per III, Fig. 2). Source and background spectra from both epochs were fit joir itly, us- 



ing M CMC simulations as implemented in the X-ray software package XSPEC (lArnaud 



19961 ). We assumed fiat priors for the absorbing column density Nh and the power law 
spectral index F, and a logarithmic prior for (absorbed) 2007 and 2009 X-ray fluxes in 
the energy range from 1 to 7 keV. The fitted absorption A^^ is 7.23(7.04,7.42) x 10^^ 
cm-2, F = 2.40(2.33,2.46), Fi_7 kev(2007) = 2.70(2.64,2.76) x lO^^^ g^gs cm'^ s-\ and 
Fi_7 kev(2009) = 2.81(2.78, 2.84) x 10^^^ ergs cm'^ g-^ (We neglect dust scattering, so the 
value of Nh here includes contributions from both absorption and scattering.) The best 
estimate of the flux increase is 4.1%, but uncertainties are large (the 90% confidence interval 
is (1.7%-6.6%), using the 0.05 and 0.95 quantiles of the MCMC draws). Since the thermal 
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contribution to the integrated spectrum is small, the increase is evidently in the synchrotron 
emission. 

The rate of X-ray flux increase is 1.7 ± 1.0% yr^^. This is consistent with our image- 
based result, and comparable to the rate of increase in the radio flux density, estimated at 
2% yr-^ in Paper II and 1.221°;^^% yr'^ bv iMurphv et all toO^ {la errors). 



Discussion 



For a distance of 8.5 kpc and a mean shock radius of about 50" (the radius of the bright 
ring in Fig. 1), our measured expansion rate gives a shock velocity of 13,000 km s~^; the 
SE-NW extensions beyond the bright ring are at about 60" for a shock velocity of almost 
16,000 km s~^. These values bracket the spectroscopically deduced velocities of order 14,000 
km s~^ (Paper IV) and indicate that our distance estimate is a good approximation. 

While our measured expansion rate gives an upper limit to the age of Gl.9+0.3 of 156 yr, 
it is likely considerably younger because of blast-wave deceleration. Significant deceleration 
is expected for either a Type la or core-collapse (CC) SN explosion. The simplest Type la 
SNR model consists of ejecta with an ex ponential ejecta density profil e, pe oc exp(— f /t>e)t"^. 



expanding into a uniform ambient ISM (iDwarkadas fc Chevalieilll998h . The velocity scale Ve 
is equal to 24:40EI(^ [Me/ Mch)~^^'^ km s~"^, where £^51 and Me/Mch are ejecta kinetic energy 
(in 10^^ ergs) and mass (in units of the Chandrasekhar mass Men)- The free expansion 
velocity of ejecta at the reverse shock is very high, V = 19600 {Rr/2 pc) (t/100 yr)"' km s-\ 
or v/ve = 8.03 {Rr/2 pc) {t/ 100 yi)^^ E^^^'^ {Me/ Mch)^^'^ , where R,. and t are the reverse 
shock radius and age of Gl.9+0.3. At very early times, the separati on between forward and 
revers e shocks is a small fraction of their radii, ~ 10% in the model of IDwarkadas fc Chevalier 
(119981 ) with the standard postshock compression r of 4 (corresponding to the adiabatic index 
7 = 5/3), and less for a cos mic ray modified blast wave with a larger compression. Using mass 
conservation ( McKeelll974j ). this fractional separation ca n be estimated as ^ (rj (r— 1))^/^ — 1. 



There is growing evidence for r > 4 in young SNR; e.g., IWilliams et al.l (1201 ll ) found a large 
(r ~ 12) blast wave compression in the LMC SNR 0519—69.0. In their hydrodynamical model 
with 7 = 1.18, the shocked region thickness is only ~ 3% at early times. We then used a 
thin-shell ap proximation (Rr = R); for ambient medium mass M varying as r^, eqs. (37) 



and (38) in IChevahen (120051 ) reduce to {[(3 + 1]?/ — 1 — 4x — 6x^) {x^dy / dx + [1 — x\y) = 
(3y {1 + 3x[l + 2x + 2x^]), where y = -^^^ exp(f/t>e) and x = ^. At very early times, 
y ^ 1 for well-behaved, asymptotic thin-shell solutions. This imposes an initial condition 
?/ = 1 at a; = 0; it is straightforward to verify that y = 1 + 4x + 6x^ is an exact solution 
of the ordinary differential equation for y that satisfies this initial condition. For uniform 
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ambient medium (/3 = 3), we arrive at an analytic, parametric solution for the dimensionless 
blast wave radius r' = R/ R' = (1 + ^Vf,/v + 6(fe/f )^)^''^ (2t>/t>e)^^^ exp(— i;/3fe), where R' = 
2.19 {Me/MchY^^ ^Q^^^ pc and uq is preshock H density in cm~^. A good match to the 
measured expansion rate of 0.64% yr~^ and the remnant's radius of 2 pc is provided by a 
fiducial SN model with E^i = 1 and Me = Mch at an age of 110 yr and with uq = 0.022 cm~^. 
For this model, v/ve = 7.3, r' = 0.25, the dimensionless time t' = r' (v /3^^'^2ve) = 0.12, 
and the deceleration parameter m = d\nR/d\nt = 0.69. 

In the CC scenario, we favor explosions of stripped compact cores over red supergiant 
explosions, because of the very high shock velocity and the presence of shocked Fe and other 
heavy elements in the outer ejecta layers. The outermost ejecta in explosions of stripped 
massive cores ( WR stars) are well descri bed by a power law density profile, pe oc v~'^t~^ 
with n = 10.2 (IMatzner &: McKed Il999l ). The blast wave radius r i nitially increase s as 
_ ^(n-3)/(n-2) g g]\j- gjgQ^g^ g^j-Q exp a nding into a WR stellar wind ( Chevalien Il982l ) . so 
that m = 0.88. IChevalier fc FranssonI ( l2006l ) find that m drops slightly to 0.866 ~ 20 yr 
after the explosion (for a representative WR progenitor), corresponding to n = 9.5 (the ejecta 
density profile becomes less steep as v decreases). Since X-ray and radio emission decreases 
with time while the blast wave propagates through the wind, Gl.9+0.3 cannot be at this 
stage of evolution. Instead, the blast wave must have passed the WR termination shock and 
now be propagating into the shocked wind with an approximately constant density. The WR 
termination shock can be less than 2 pc cl ose to the Galac t ic cen ter where the ISM pressure 
is much higher than in the solar vicinity. IChevalier et al.l ( 120041 ) investigated the evolution 
of bubbles blown by WR stars in high-pressure environments, and found that the wind 
termination shock Rt stalls at Rt ~ l.Svl^'^M^^^ {p/k/10^ cm^^K) pc for ISM pressures 
p/k higher than 2.5 x 10^M_5t>8 {t/twR)~'^ cm^^ K (fg, M-s, and t are the WR wind speed, 
mass loss rate, and duration in units of 1000 km s^^, 10~^Mq yr~^, and ti^R = 3 x lO'* 
yr, respectively). If the current blast wave radius is much larger than the wind termination 
shock, then m = {n — 3)/n (or 0.68 for n = 9.5; Chevalier 1982), comparable to our estimates 
for the Type la scenario. Intermediate values between 0.68 to 0.87 arise if the termination 
shock radius is a sizable fraction of the current remnant radius. This will affect the estimated 
age. With m = 0.7 (0.8), Gl.9+0.3 is 110 (125) yr old. Ksenofontov et al. (2010) modeled 
Gl.9+0.3 with power-law ejecta with n = 7, obtaining an age of only 80 yr (they also relied 
on estimates of the remnant's expansion from Papers I and II). 

Reahstic high-velocity ejecta density profiles in Type la SNe may deviate substantially 
from the idealized exponential profile. High velocity features are frequently seen in Type la 
SN spectra at early times, mostly in Ca II absorption, suggesting the presence of density 
structures of unknown origin. Type la explosion models differ significantly in predicted 
density profiles at high velocities, including variations in steepness and the presence or 
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absence of distinct densit y structures. T 



lis st rongly affects the deceleration of the blast 



wave, as demonstrated by iBadenes et al.l (l2003l . Figs. 2 and 3e). It is then highly desirable 



to measure the deceleration directly with present and future Chandra observations. This may 
be feasible with observations spanning the anticipated Chandra hfetime (we hope at least 
one more decade). Assuming no deceleration, the current expansion age of 156 yr at the 
mean epoch 2008 would increase to 166 yr at epoch 2018; alternatively, the expansion rate 
would decrease to 0.603% yr ~^ from the current value of 0.642% yr~^. With m = 0.7, we 
expect a smaller rate of 0.588% yr~^ in 2018. It may be possible to measure this difference, 
since significant improvements are expected with upcoming long Chandra observations, on 
account of much better photon statistics and the increased time baseline. 

The flux increase we report of about 1.7% yr^^ is consistent with both our estimate 
of the rate of increase of radio flux (about 2%; Paper II), and that reported by Murphy et 
al. (2008), 1.22+°:?^% yr~^ As discussed in Paper II, radio flux increases are difficult but not 
impossible to produce; in the absence of increasing efficiency of magnetic-field amplification 
and/or particle acceleration, they typically require expansion into a uniform medium. We 
should note that prompt radio emission from supernov ae is observed to dr op rapidly with 



time at frequencies at which it is optically thin (e.g., iPanagia et al.l 120061 ): this emission 
is typically modeled as synchrotron emission from a region occupied by magnetic field and 
relativistic electrons, with a constant fraction of shock energy going into each component, 
as a shock moves through a dense stellar wind. The increasing radio and X-ray flux from 
Gl.9+0.3 argues against current interaction with a stellar wind. In the formalism of Katsuda 
et al. (2010), Appendix A, for expansion into a uniform medium, the intensity at radio 
wavelengths /j, varies at the fractional rate I~^dly/dt = p/t, wherep = m(4+a) — (3 + a) and 
Jj, oc i>~°'. In that work, the relativistic-electron density and magnetic- fleld energy density 
are both assumed proportional to the post-shock energy density pM^h- Then for a = 0.62, 
appropriate for the integrated radio spectrum of Gl.9+0.2, we have p = 4.62m — 3.62. The 
integrated flux S^, oc luR^, giving a fractional rate 71 = S'^dS^/dt = {p + 2m) /t. For 
Gl. 9-1-0. 3, we consider two cases: undecelerated expansion (m = 1 and t ~ 156 yr), and 
m = 0.7 (t ~ 110 yr). The undecelerated case then gives a fractional flux increase rate of 
about 2% yr^^, while the decelerated case gives 1% yr~^, so that the radio results nicely 
bracket the theoretical expectations. 

For a power-law electron distribution with an exponential cutoff above some -Emax, the 
rate of change of intensity and flux will be different in general for observing frequencies 
V > Vf. = 1.82 X IQ^^E'^^^B, the characteristic synchrotron frequency emitted by electrons 
with energy -Emax in a magnetic fleld B (e.g., Pacholczyk 1970; cgs units). If rises with 
time, the X-ray flux will rise faster (or decline more slowly) than the radio, while if Vc drops 
with time, the reverse is true. So the question of whether the X-ray rate of change is larger 
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or smaller than the radio is of considerable theoretical interest. The same Appendix from 
Katsuda ct al. (2010) gives the additional fractional rate of change of intensity (or flux) as 
AT?. = {a/2t) yjv jvc-, where oc t". For an average over the spectrum of Gl.9+0.3, we take 
vjvc = 1.5 (or hu = 3.3 keV, about the peak in the integrated spectrum). The behavior 
with time of Uc depends on the mechanism limiting electron acceleration. If radiative losses 
limit it, Katsuda et al. (2010) show that a = 2m — 2, while if the finite remnant age is 
the limitation, a — 7{m — 1) + 2. Then undecelerated motion predicts the same rate of 
change in X-rays as radio for loss-limited acceleration, and a faster rate, about 2.8% yr~^, 
for age-hmited acceleration. For decelerated expansion, both possibilities give slower rates 
of increase: about 0.6% yr~^ for loss-limited, but only slightly less than 1% yr"-*^ for age- 
limited acceleration. Subsequent observations should allow both the radio and X-ray rates 
of flux increase to be determined more accurately, offering the possibility of determining 
the mechanisms limiting particle acceleration in Gl.9+0.3. We point out that age-limited 
acceleration would have to limit ion acceleration as well, while if electrons are loss-limited, 
ions might well be accelerated to much higher energies than those we deduce from electron 
synchrotron emission. 

6. Conclusions 

1. The change in morphology of the X-ray emission from Gl. 9-1-0. 3 is well- described by 
uniform expansion at a rate of 0.642 ± 0.049% yr~^ between 2007 and 2009, consistent 
with earlier observations. Undecelerated expansion would imply an age of 156 ± 11 
yr, but reasonable expectations for either a Type la or core-collapse event suggest 
deceleration with R oc giving an age of about 110 yr or an explosion date of 1900. 

2. The total X-ray flux is increasing at the rate of 1.7±1.0% yr^^ (90% confldence limits), 
a rate comparable to radio rates. This increase makes it highly unhkely that Gl. 9-1-0. 3 
is expanding into a stellar wind. 

3. Direct determination of deceleration may be possible within the lifetime of Chandra. 
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Table 1. Best-Fit Parameters 



Parameter 


Value 


Expansion Rate 


(0.642±0.049)% yr-i 


Surface-Brightness Factor 


1.01 ±0.02 


Flux Increase Rate (image) 


(1.89±0.86)% yr-i 


Flux Increase Rate (spectrum) 


(1.7 ± 1.0)% yr-i 


Expansion Center RA 


17:48:45.570±0.005 


(J2000) Dec 


-27:10:06. 94±0.08 



Note. — Uncertainties arc 90% confidence intervals. Errors 
on the expansion center are statistical only; Chandra pointing 
uncertainties are greater. 
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R.A. a2000) 



Fig. 1. — Chandra X-ray image from 2009, smoothed with platelets ( IWillettI 120071 ). The 
fitted expansion center is indicated by the + sign. Surface brightness is in cts px~^. 
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Fig. 2. — Profiles from 2007 (red) and 2009 (blue) along position shown in Figure 3. Hori- 
zontal scales are distance from the expansion center in arcseconds. Top: SE limb; bottom, 
NW limb. Expansion is evident in each. 
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Fig. 3 
taken. 



— Raw 2009 Chandra image, showing region from which the profiles of Fig. 2 were 



